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ABSTRACT 

Context. It has recently been suggested that the compilation of the calibrated time-ordered-data (TOD) of the Wilkinson 
Microwave Anisotropy Probe (WMAP) observations of the cosmic microwave background (CMB) into full-year or multi- 
year maps may have been carried out with a small timing interpolation error. A large fraction of the previously estimated 
WMAP CMB quadrupole signal would be an artefact of incorrect Doppler dipole subtraction if this hypothesis were 
correct. 

Aims. Since observations of bright foreground objects constitute part of the TOD, these can be used to test the 
hypothesis. 

Methods. Scans of an object in different directions should be shifted by the would-be timing error, causing a blurring 
effect. For each of several different timing offsets, three half-years of the calibrated, filtered WMAP TOD are compiled 
individually for the four W band differencing assemblies (DA's), with no masking of bright objects, giving 12 maps 
for each timing offset. Percentiles of the temperature-fluctuation distribution in each map at HEALPix resolution 
A^side = 2048 are used to determine the dependence of all-sky image sharpness on the timing offset. The Q and V bands 
are also considered. 

Results. In the W band, which is the band with the shortest exposure times, the 99.999% percentile, i.e. the temperature 
fluctuation in the ~ 503-rd brightest pixel, is the least noisy percentile as a function of timing offset. Using this statistic, 
the hypothesis that a —25.6 ms offset relative to the timing adopted by the WMAP collaboration gives a focus at least as 
sharp as the uncorrected timing is rejected at 4.6a significance, assuming Gaussian errors and statistical independence 
between the maps of the 12 DA/observing period combinations. The Q and V band maps also reject the —25.6 ms 
offset hypothesis at high statistical significance. 

Conclusions. The requirement that the correct choice of timing offset must maximise image sharpness implies that 
the hypothesis of a timing error in the WMAP collaboration's compilation of the WMAP calibrated, filtered TOD is 
rejected at high statistical significance in each of the Q, V and W wavebands. However, the hypothesis that a timing 
error was applied during calibration of the raw TOD, leading to a dipole-induced difference signal, is not excluded by 
this method. 

Key words, cosmology: observations - cosmic background radiation 



1. Introduction 

■ The nature of the large-scale cosmic microwave background 
, (CMB) signal in the Wilkinson Microwave Anisotropy 

■ Probe (WMAP) observations (Bennett et al. 2003b) is of 
fundamental importance to observational cosmology. A lack 

. of structure on the largest scales in maps of CMB temper- 
' ature fluctuations would be a hint that the comoving spa- 
' tial section of the Universe is multiply connected (de Sitter 
1917; Friedmann 1923, 1924; Lemaitre 1931; Robertson 
1935) and that its size may have been detected (Starobinsky 
1993; Stevens et al. 1993). The WMAP observations have 
confirmed that the large-scale CMB signal analysed either 
as a quadrupole signal or as a two-point autocorrelation 
function signal is weak (Spergel et al. 2003; Copi et al. 2007, 
2009). Several analyses indicate that either a comoving 
space (Spergel et al. 2003; Aurich et al. 2007; Aurich 2008; 
Aurich et al. 2008, 2010) or a Poincare dodecahedral space 
(S^/r) model (Luminet et al. 2003; Aurich et al. 2005a, 
2005b; Gundermann 2005; Caillerie et al. 2007; Roukema 
et al. 2008a, 2008b) fit the WMAP data better than an 
infinite, simply connected flat model, while other analyses 



disagree (Key et al. 2007; Niarchou & Jaffe 2007). Highly 
significant evidence for either an infinite or a finite model 
of comoving space has not yet been obtained. 

Liu & Li (Liu & Li 2008; Li et al. 2009; Liu & Li 2009; 
Liu et al. 2010) have suggested that several systematic er- 
rors may have been made in the data analysis pipelines 
used by the WMAP collaboration in the compilation of 
their original time-ordered-data (TOD) into single-year or 
multi-year maps. Liu & Li (2008) recommended that new 
maps should be made using the full TOD in order to avoid 
some of these suspected errors. This requires moderately 
heavy computational resources (RAM and CPU) on typ- 
ical present-day desktop computers. Aurich et al. (2010) 
sidestepped this by using a phenomenological correction for 
the hot pixel effect (Liu & Li 2008). They found that their 
correction applied to the 5-year WMAP W band data re- 
moves the anti-correlation in temperature fiuctuations at 
angular separations of nearly 180° and improves the fit of 
a model to the data. 

Liu & Li (2009) carried out a full analysis of the 5-year 
WMAP TOD. Their analysis pipeline appeared to give test 
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Fig. 1. Dipole-induced difTerence all-sky map using the WMAP quaternion history of the first year of WMAP observations and 
5t = in the V band, i.e. a timing offset of —38.4 ms (see Sect. \'2.'2} . shown in the hybrid, cylindrical equal-area (Lambert 1772), 
interrupted CoUignon (Collignon 1865) spherical projection, i.e. the HEALPix {H = 4, K — 3) spherical projection (Roukema & 
Lew 2004; Calabretta & Roukema 2007; Gorski et al. 2005), centred on the Galactic Centre, with latitude increasing upwards and 
longitude increasing to the left. The grey scale is linear from — 10/iK (black) to +10fiK (white). No sky observations are used, i.e. 
only Liu et al. (2010) 's Eqs (1) and (2) as encoded in their software are used to process the quaternion data. Standard sky masks 
are used. This is the 80-th map iteration. 



results largely compatible with those of the WMAP collabo- 
ration. However, the CMB quadrupole amplitude was found 
to be much weaker, and sub-degree power was also found 
to be a little weaker. Later, Liu et al. (2010) traced the 
difference between their analysis and that of the WMAP 
collaboration to a difference in interpolating the timing 
of individual observations from the times recorded in the 
TOD files. The TOD FITS format files contain observa- 
tional starting times in both the "Meta Data Table" and 
the "Science Data Table". The full set of observing times 
of individual observations are not recorded in either data 
table; they must be interpolated from the smaller set of 
observing times that are recorded. The starting times in 
the Meta Data Table are 25.6 ms greater than the cor- 
responding times in the Science Data Table. A computer 
script for checking this is given in Appendix [Xj The dura- 
tions of individual observations in the Q, V and W bands 
are 102.4 ms, 76.8 ms, and 51.2 ms, respectively (Section 
3.2, Limon et al. 2010). The 21 April 2010 version of the 
WMAP Explanatory Supplement (Section 3.1, Limon et al. 
2010) discusses the relation between the start times in the 
two tables, but does not refer to the 25.6 ms offset, i.e. half 
of a W band observing interval or a quarter of a Q band 
observing interval. Bennett et al. (2003a) state (Sect. 6.1.2) 
that "a relative accuracy of 1.7 ms can be achieved between 
the star tracker(s), gyro and the instrument", i.e. the com- 
bined random and systematic error should be much smaller 
than 25.6 ms. Given that the timing offset is numerically 
implicit in the WMAP TOD files, it could, strictly speak- 
ing, be referred to as the WMAP collaboration's "implicitly 
claimed but ignored" timing offset. However, for simplicity, 
the timing offset will be referred to here as Liu et al. (2010) 's 
hypothesis. 



Using their Eqs (1) and (2), Liu et al. (2010) showed 
that the error in estimating the Doppler dipole of the space- 
craft velocity implied by the would-be timing error can 
be used to produce a map (Fig. 2 left, Liu et al. 2010) 
that matches the direction and amplitude of the CMB 
quadrupole signal as estimated by the WMAP collabora- 
tion (Fig. 2 right, Liu et al. 2010) remarkably well, de- 
spite using only the history of spacecraft attitude quater- 
nions (representing three-dimensional directions) and no 
sky observations. This is easy to check using the WMAP 
calibrated, filtered, 3-year TOD and Liu et al. (2010) 's 
software. Figure [1] shows a dipole-induced difference mapQ 
which is consistent with Liu et al. (2010) 's dipole-induced 
difference map shown in their Fig. 2 (left), for a half- 
observing-interval offset in the V band, i.e. —38.4 ms. 
Using a velocity dipole model and a —25.6 ms offset, Moss 
et al. (2010) find a similar but weaker effect of reduc- 
tion in the quadrupole amplitude. Earlier, Jarosik et al. 
(2007) (Sect. 2.4.1) noted that an error during the calibra- 
tion process could introduce an artificial, dipole-induced 
quadrupole. 

Interpolation of the observing times and sky positions 
from the WMAP TOD files cannot be avoided. Given that 
the WMAP Explanatory Supplement does not presently 
account for the 25.6 ms offset between the two data tables 
in each TOD file, and given the potential importance of a 
70-80% overestimate of the CMB quadrupole amplitude, 
it is useful to see if a relatively simple, robust method of 



^ Liu et al. (2010) refer to this as a "differential dipole" map. 
However, this term is frequently used to refer to monopole- 
subtracted dipole maps. Moreover, the difference between two 
(fixed) dipole signals in different directions is itself also a dipole, 
not a quadrupole. To avoid misinterpretations, the term "dipole- 
induced difference" is suggested here. 
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Fig. 2. W band cylindrical equal-area projection of the 
Galactic Centre using 199 days of the 2002 WMAP TOD 
(Sect. 12. ip . compiled using a minor modification of Liu et al. 
(2010)'s analysis script (Sect. [2]2}, for a fractional timing off- 
set St = —5 (i.e. an offset of —256 ms relative to that of Liu 
et al.). Galactic latitude increases upwards, galactic longitude 
increases to the left. The image size is 5.8° x 1.7°, centred at the 
Galactic Centre. Individual pixels are 1.7' in side length. The 
temperature-fluctuation grey scale ranges from black (-20 mK) 
to white (-1-40 mK). 



analysing the TOD can determine which of the two timing 
interpolation methods, that of the WMAP collaboration or 
that proposed by Liu et al. (2010), is correct. 

The WMAP TOD include many observations of bright 
foreground objects, including Solar System planets and ob- 
jects in the Galactic Plane. The complex scanning pattern 
of WMAP implies that individual objects are likely to be 
scanned in different directions over a series of successive ob- 
servations. The greatest differences in scanning directions 
are likely to occur over observing intervals of weeks or more. 
An error in the assumed sky direction due to a timing er- 
ror should cause observations of a compact source made 
at different times to be slightly offset from one another in 
the compiled map, i.e. it should cause an overall shift in sky 
position and a blurring effect. For a large enough timing off- 
set, double imaging of compact sources should occur (pub- 
lic communication, Crawford 201CQ; see also Moss et al. 
2010). The W band observations, with the highest angular 
resolution, have an estimated FWHM beam size of about 
12-12.6' (Table 5, Page et al. 2003). An error of 25.6 ms 
corresponds^ to an error of about 4.0'. Hence, a slight blur- 
ring of the images of compact sources rather than double 
imaging should occur. If the WMAP collaboration's tim- 
ing offset were correct, then this blurring could explain the 
slight loss of sub-degree power found by Liu et al. (2010) as 
an artefact for the latter's choice of timing offset. The blur- 
ring could be detectable in beam shape analysis, but this 
may not be easy. Page et al. (2003) model the beam shape 
based on the physical structure of the WMAP receivers, but 
estimate centroids iteratively. It is not obvious that this it- 
erative procedure would have detected an artificial blurring 
effect of just a few arcminutes. 

Nevertheless, the blurring effect should be statistically 
detectable in maps compiled from the TOD. The effect 
should be strong if the timing error is exaggerated beyond 
that proposed by Liu et al. (2010). By considering the tim- 



^ http : //cosmocof f ee . inf o/viewtopic.php?t=1537 
^ One spin during 2.2 minutes, i.e. 129.3 s, (e.g. Sect 3.4.1, 
Hinshaw et al. 2003) corresponds to 360° sin(70.5°) — 339° in 
great circle degrees, so 25.6 ms corresponds to 4.0'. 



Fig. 3. Image of the Galactic Centre, as in Fig. [J] for a frac- 
tional timing offset St = (proposed by Liu et al. 2010). 




Fig. 4. Image of the Galactic Centre, as per Fig. [2l for a frac- 
tional timing offset St = 0.5 (WMAP collaboration value). 



*■ ... * 



Fig. 5. Image of the Galactic Centre, as per Fig. [2l for a frac- 
tional timing offset St = 5. 



ing offset as a free parameter, it should be possible to find 
the value that minimises blurring. This should correspond 
to the correct choice of timing interpolation. In principle, 
since a timing error also induces an error in absolute sky 
position, this could also be used to determine the correct 
method. However, it is conceivable that a positional offset 
may have been at least partially removed at some stage in 
the pointing calibration, making it hard to detect, as sug- 
gested by Moss et al. (2010). On the other hand, a blurring 
effect is unlikely to have been removed at any step in the 
production of WMAP maps. Hence, the analysis here uses 
the blurring effect. In Sect. [H an empirical method of esti- 
mating the blurring using all-sky maps is presented. Images 
of some Galactic Centre objects as a function of timing off- 
set and the results of the statistical analysis are presented 
in Sect. |21 Conclusions are given in Sect. HI 
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2. Method 

2.1. Observational data files and differencing assemblies 

Compilation of TOD into maps requires a moderately heavy 
usage of RAM, CPU and disk space on present-day desktop 
computers. For this reason, only a subset of the available 
WMAP TOD is used rather than the full TOD. From each 
of the years 2002, 2003 and 2004 the files for the first (by 
filename) 199, 199, and 198 daysQ respectively, of filtered, 
calibrated TOE0 are used, i.e. slightly over six months of 
data during each year. This allows a convenient filename 
pattern. Liu & Li (2009) analysed the unfiltered, calibrated, 
5-year WMAP data release. The filtering should not signif- 
icantly affect the questions of interest here. 

The primary analysis is done here using the W band, 
since this has the best resolution. There are four W differ- 
encing assemblies (DA's), giving 12 TOD sets that in many 
ways can be considered to constitute independent observa- 
tions. Complementary analyses using the Q and V DA's, 
with 6 TOD sets each, are also carried out. 

2.2. Compilation of the TOD 

Liu et al. (2010) have made their script for compiling the 
TOD publicly available^ For the purposes of the present 
analysis, minor modifications of this script and associated 
libraries have been made in order to run the script using the 
GNU Data Language (GDL)13 In particular, the data flag 
parameter daf_mask has been set to 1 in order to include 
all planet observations when compiling the TOD, and the 
processing mask has been set to in order not to hide any 
bright objects from the analysis. 

Moreover, the center parameter for determing the tim- 
ing offset is generalised to an arbitrary floating point value. 
Here, this is written St, i.e. the timing offset expressed as a 
fraction of an observing time interval in a given band. This 
is used as a fraction for interpolation between quaternions. 
Thus, (5t = is Liu et al. (2010)'s primary hypothesis of 
what is the correct offset, and St = 0.5 is what the WMAP 
collaboration believes to be correct. In the W band, St — 
corresponds to —25.6 ms relative to St — 0.5. In order to 
reduce computational time, map iteration is not done. 

The planets are useful for this analysis because of their 
high brightness. However, some of them (especially Jupiter) 
move by up to a few degrees during the times when they 
happen to be scanned. In directions orthogonal to their 
movement, a blurring effect is to be expected for a wrong 
choice of timing offset. In directions parallel to their move- 
ment, a wrong choice of timing offset can either cause blur- 
ring, or cause some observations that should be mapped to 
different pixels to incorrectly be mapped to the same pixel, 
i.e. creating a false focussing effect. The degree to which the 
latter effect is important depends on the degree to which 

* The year 2004 has only 198 days because the file that would 
have the label "20041482358_20041492358" is absent from the 
data set. 

^ |http : //lamb da .gsfc. nasa.gov/product/map/dr2/ I 
|tod_f cal_get . cfi? 

http : //cosmoc of f ee . inf o/viewtopic.php?p=4525[ 

http : //dpc . aire . org. cn/data/wmap/09072731/release_vT7] 
LSOurce_code/vl/ 

^ |http : //co smo . torun . pl/GPLdownload^^ 

|LLmapmaking_GDLpatches/LLmapmaking_GDLpatches_0 .0.3. tbz| 



the observations of a given planet are discontinuous and 
the way that this discretisation interacts with the scanning 
pattern directions. In practice, semi-quantitative examina- 
tion of images of Jupiter suggests that the false focussing 
effect is present, but weak, for some values of St. An anal- 
ysis using Jupiter alone would need to take this effect into 
account. 

2.3. Pixel resolution and the blurring effect 

Since the blurring effect expected from a 25.6 ms timing 
error would be on a scale of ~ 4.0', the commonly used 
HEALPix resolution (Gorski et al. 1999) of N^Me = 512 
would be insufficient, since this gives pixel sizes of 6.9'. 
However, very high resolutions would give very large file 
sizes. As a compromise, iVgido = 2048 is adopted here, giv- 
ing a pixel size of 1.7' and one-dimensional FITS binary 
table file sizes of about 190 Mib, or two-dimensional FITS 
image file sizes of about 400 Mib, using the the HEALPix 
{H = A, K = 3) spherical projection (Roukema & Lew 
2004; Calabretta & Roukema 2007; Gorski et al. 2005). The 
utility program HPXcvT from the WCSLIB library can be 
used to convert from one-dimensional to two-dimensional 
formats Image files that use the HEALPix spherical pro- 
jection can be viewed using standard astronomical FITS 
image viewing tools for visual checking of the blurring ef- 
fect. 

Figures [5H5] show the Galactic Centre from the 2002 
WMAP TOD compiled as discussed above, for St — 
—5,0,0.5, and 5. Figures [3] and |4] show the timing options 
of Liu et al. (2010) and the WMAP collaboration, respec- 
tively. The blurring effect is not obvious by visual inspec- 
tion. These two images appear about equally sharp. In con- 
trast. Figs [5] and [HI in which the timing offsets are greater, 
clearly show both a positional offset and blurring. 

2.4. Temperature-fluctuation distribution percentiles 

The blurring in Figs [2] and [5] is the key to the method used 
here. The numbers of very bright pixels in these two figures 
are much lower than in Figs [3] and HI This suggests a statis- 
tic to characterise the degree of sharpness of the image: the 
number of pixels N\i{ST /T,St) above a given temperature 
fluctuation threshold ST /T in a map calculated at a given 
timing offset St. Maximising iVb over St should give the map 
most likely to be correct. 

The brighter the threshold ST/T, the fewer the num- 
ber of pixels that will lie above this threshold, increasing 
Poisson noise. Thus, fainter thresholds ST /T should give 
stronger signals. 

The situation for intermediate brightness pixels is less 
simple. The signal in the very bright pixels in Figs |3] and |4] 
has been spread out into intermediate brightness pixels in 
Figs m and increasing the numbers of the latter. However, 
the signal in regions of intermediate brightness in Figs [3] 
and 0] has also been blurred out of the same brightness 
interval, down to fainter levels. For typical profiles of as- 
tronomical objects, the number of pixels in each different 
brightness range should increase as the threshold decreases 
while remaining well above zero. Hence, the loss of pixels 
that become blurred below a given threshold ST jT should 
outweigh the gain in pixels from intrinsically higher bright- 

* |http : //www ■ atnf . csiro . au/people/mcalabre/WCS/l 
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Table 1. Percentile noisiness: mean of as 
band. 



(over St) in the W 



p 99% 99.9% 99.99% 99.999% 99.9999% 
{crs^)st 0.056 0.041 0.038 0.023 0.051 




Fig. 6. Mean normalised temperature-fluctuation percentiles 
Sp = (sp) [Eq. ©I in compilations of W band WMAP TOD, 
for pixel fractions p = 99.99%, 99.999%, and 99.9999%, as a 
function of fractional timing offset St. The values St — and 
St = 0.5, favoured by Liu et al. (2010) and the WMAP collab- 
oration, respectively, are shown as vertical lines. The error bars 
are standard errors in the mean over the 12 DA/year com- 
binations. For clarity, offsets in St of ±0.05 are applied to the 
p = 99.99% {+) and' 99.9999% (x) data points. 



nesses. That is, in general, an incorrect value of the frac- 
tional timing offset St should decrease at any positive 
value of 6T/T > 0. 

Approaching low (positive) temperature-fluctuation lev- 
els, the measurement noise in ST/T should proportion- 
ally increase, so that noise in A^b increases. The effects 
of smoothing with negative fluctuations will also arise at 
low positive temperature- fluctuation levels, reducing the ef- 
fectiveness of testing TVb. Hence, the least noisy threshold 
6T/T for optimising the sharpness of images in the map 
should lie somewhere between the extremes of low positive 
ST/T and high 6T/T. 

Rather than choosing fixed values of 6T/T and study- 
ing Nt{ST/T,St), the approach adopted here is to con- 
sider ST/T{pN, 6t), where p is a fraction close to unity and 
N = 12(2048)^ is the number of pixels. That is, the per- 
centiles of the (sorted) temperature-fluctuation distribution 
are determined for several values of p, as a function of the 
fractional timing offset St. 

Use of different DA's and different (half-)years of data 
should give estimates of the random error at any given p and 
St. Averaging over 6t at a given p makes it possible to find 
the fraction p that gives the least noisy percentile. Thus, at 
a given DA and year combination i and a percentage p ^1, 
the percentile is a function 



(1) 



Table 2. Mean difference in temperature-fluctuation per- 
centiles A{St') ~ ( {ST/T)p{St = 0.5) - {ST/T)p{St') ) for 6 
DA/year combinations in each of the Q and V bands in mK, 
where p = 99.999%, and standard error in the mean and signif- 
icance level.'' 



band 


A{St' = 


0) A(0.167'') 


A(0.25") 




A/CTA 


Q 


1.44 






0.05 


29 


Q 






0.44 


0.08 


5.5 


V 


0.70 






0.06 


11.7 


V 




0.28 




0.06 


4.8 



''The time offsets St and St' are in units of one observing 
interval, i.e. 102.5 ms or 76.8 ms for Q or V, respectively. 
''A shift of -25.6 ms relative to the WMAP collaboration 
choice in the V band. 

'^A shift of —25.6 ms relative to the WMAP collaboration choice 
in the Q band. 



within each sample i at each p. That is, {5T/T)p is evalu- 
ated on a series of maps calculated for — 5 < <5t < 5, attain- 
ing a minimum and a maximum bp. This is normalised 
to 



sliSt) 



(ST/T); (St) - a; 



(2) 



In order to give approximately equal weight to the different 
samples i, the percentiles are normalised over St internally 



The mean and standard error in the mean of Sp over the 
DA/year combinations i can then be estimated at each St. 
These are written Sp and , respectively. 



3. Results 

3.1. W band 

Calculations on 4-core, 2.4 GHz, 64-bit processors 
with 4 Gib RAM, using GDL-0.9~rc4 running under 
GNU/Linux, took about 3 hours per map. Table [T] lists 
the mean of tTsp (St) at several values of p for the W band 
analysis. This shows that p — 99.999%, i.e. the fraction 
p which selects the « 503-rd brightest pixel in any given 
map, gives the least noisy percentile. Figure [5] shows the 
mean normalised temperature-fluctuation percentiles Sp for 
this and two other values of p. The 599.999% threshold 
shows a sharp maximum at St = 0.5. The higher per- 
centile, i.e. at p — 99.9999%, is noisier, but still signifi- 
cantly favours St = 0.5 over any other value of St evalu- 
ated. The p = 99.99% percentile is also noisier than that 
at p = 99.999%, and marginally favours St = 0.75 over 
St = 0.5. Liu et al. (2010)'s hypothesis oi St = is not 
favoured at any fraction p. 

Exact symmetry around the correct value of St is not 
necessarily expected, since the scan paths are curved (non- 
geodesic) paths on the 2-sphere, not straight lines in the 
2-plane, and the source objects are not isolated from one an- 
other. However, approximate symmetry can reasonably be 
expected. It is obvious that Fig. [5] is not symmetric around 
St ^ 0. Hiding the part of the figure to the left oi St = — 1 
shows that the figure does appear to be approximately sym- 
metric around St = 0.5. 

A quantitative estimate of a lower limit to the signifi- 
cance to which the St = hypothesis can be rejected, given 
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that it must maximise Sp{dt), can be obtained by testing 
the hypothesis that s'p{St = 0) > Sp{St = 0.5) under the 
assumption that the error distributions are Gaussian. That 
is, the hypothesis that that a —25.6 ms offset relative to 
the timing adopted by the WMAP collaboration gives a 
focus at least as sharp as that obtained using the latter's 
timing choice is tested. The percentiles shown in Fig. [5] are 
Sp(0) = 0.901 ± 0.019 and Sp(0.5) = 0.989 ± 0.005, where 
p = 99.999%, i.e. where Sp is the percentile that gives the 
least noisy result, so the hypothesis is rejected at 4.6cr sig- 
nificance. 

3.2. Q and V bands 

A small number of maps were calculated in the Q and V 
bands, which have lower angular resolution and fewer num- 
bers of DA's. Only a few maps were calculated, so the per- 
centiles are not normalised. Instead, their differences are 
calculated. The Liu et al. (2010) hypothesis may be inter- 
preted in the Q and V bands either as a fractional time 
offset St ^ (where St = 0.5 is the WMAP collabora- 
tion's choice), or as a fixed offset in time units of —25.6 ms, 
since this is the undocumented difference between the start 
times in the Meta Data Set and the Science Data Set in 
the TOD files. The latter difference should be more diffi- 
cult to test than the former, since the interval is smaller 
in both cases. These differences are calculated individually 
in each DA/year combination. Table [5] shows their statis- 
tics. The weaker form of Liu et al. (2010) 's hypothesis is 
rejected at very high significance (> 4.8(7) in all cases. In 
other words, the temperature fluctuation in the Ri503-rd 
brightest pixel is lower for Liu et al. (2010) 's hypothesised 
offset compared to its value for the WMAP collaboration's 
hypothesised offset at very high significance in all cases. 

4. Conclusion 

The difference in image sharpness for Liu et al. (2010) 's and 
the WMAP collaboration's respective timing interpolation 
choices are not obvious in the images of the Galactic Centre 
shown in Figs |3] and HI However, use of the all-sky images 
gives estimates of image sharpness that depend strongly 
on the fractional time offset St. These calculations do not 
require any fitting of source or beam profiles, nor any selec- 
tion of one or more preferred sources, nor any assumptions 
regarding pointing accuracy. The data analysis pipeline is 
performed using Liu et al. (2010) 's script, with only minor 
modifications, on slightly under a quarter of the full set of 
WMAP calibrated, filtered TOD. Liu et al. (2010)'s timing 
error hypothesis is rejected at very high significance in the 
Q, V and W wavebands separately. 

The lack of an explanation in the WMAP Explanatory 
Supplement (Section 3.1, Limon et al. 2010) for the 25.6 ms 
offset between the start times in the Meta Data Set and the 
Science Data Set in the TOD files is an unfortunate gap in 
the documentation, but it clearly did not lead to an error 
of this magnitude in the WMAP collaboration's compila- 
tion of the calibrated, filtered TOD. A possible explanation 
for the offset could be that it is a convention inserted by 
software into the Meta Data Set that was forgotten about, 
rather than a physical time offset. In principle, it would be 
good if the WMAP collaboration could find the line(s) in 
their software where the 25.6 ms offset is inserted, in order 



to confirm that its origin is fully understood, and that it 
has no unintended consequences. 

Nevertheless, the fact that Liu & Li (2009) acciden- 
tally applied a time offset that approximately cancels the 
St = 0.5 CMB quadrupole in amplitude and direction, and 
that this can be modelled as a Doppler-dipole-induced dif- 
ference signal (Eq. (2), Liu et al. 2010, "deviation of differ- 
ential dipole"), remain interesting coincidences. The rela- 
tions between the quadrupole, octupole and ecliptic (e.g., 
Copi et al. 2006; Sarkar et al. 2010) have long been known. 
What remains unestablished is whether these are coinci- 
dences or artefacts. Moss et al. (2010) discuss the possibil- 
ity of relationships between the various coincidences and 
suggest that another effect might couple with the WMAP 
spacecraft velocity dipole in order to produce an artificial 
component of the quadrupole. 

One interesting possibility for further work would be 
to consider a variation on the hypothesis proposed by Liu 
et al. (2010): could it be possible that the would-be tim- 
ing error was introduced when calibrating the uncalihrated 
TOD? Although the velocity dipole is only removed from 
the data during the mapmaking step, it is first used for cal- 
ibrating the raw data into a calibrated signal in mK. So, a 
small timing error during the calibration step could corre- 
spond to artificially adding a dipole-induced difference sig- 
nal. During the following step of compiling the calibrated 
TOD into maps, the dipole would then be subtracted at 
the correct positions, but without compensating for the in- 
correct calibration. As noted above, Jarosik et al. (2007) 
(Sect. 2.4.1) have earlier considered the possibility that a 
calibration error (not necessarily induced by a timing error) 
could introduce an artificial, dipole-induced quadrupole. 

The calibration method for the 7-year WMAP analy- 
ses are stated (Jarosik et al. 2010) to be those used for 
the 5-year analyses, i.e. as described in Sect. 4 of Hinshaw 
et al. (2009), retaining the estimate of about 0.2% abso- 
lute calibration error. Since the dipole amplitude is about 
3 mK, a 0.2% error corresponds to about 6 /j,K. The 7- 
year quadrupole estimated for St — 0.5 is a/ (3/7r)C2 ~ 
14/iK (Section 4.1.1, Jarosik et al. 2010), only a factor of 
two higher than this calibration uncertainty (though the 
model-dependent systematic error is high for the infinite 
flat model). Liu & Li (Sect. 4, Fig. 6, 2009) estimate the 
quadrupole-like difference between their Ql DA St = Q map 
and the WMAP collaboration's St = 0.5 map to have an 
r.m.s. of 6.6 /xK. Moreover, the spacecraft's velocity is as- 
sumed by Hinshaw et al. (2009) and Jarosik et al. (2010) 
to be known exactly, and fits are made over intervals "typ- 
ically between 1 and 24 hours". It is not obvious that a 
timing error causing the spacecraft's velocity vector to be 
slightly misestimated, which in turn causes a quadrupole- 
like dipole-induced difference signal over a year's observa- 
tions, of amplitude ^ 5-10/iK, would have been detected 
or removed during the calibration process. Indeed, Cover 
(2009) claims the presence of considerable systematic un- 
certainty in the calibration of the WMAP TOD. 

The calibration only applies to (differential) radio flux 
density estimates, not to positional data, so it would not af- 
fect the latter, i.e. the spacecraft attitude quaternions. The 
Meta Data Tables containing the quaternions appear to be 
identical between the "uncalibrated" and "calibrated, fil- 
tered" WMAP 3-year TOD files, as can be expected. Hence, 
an incorrect calibration that leads to a dipole-induced dif- 
ference signal would not cause an arcminute-scale blurring 
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effect nor a positional error in point sources. It would not be 
detectable by the method used in this paper. Since a large 
fraction of the brightest point sources at WMAP frequen- 
cies are variable, it may not be easy to use point sources to 
calibrate the uncalibrated data independently of the dipole 
calibration. Liu & Li (2009) and Liu et al. (2010) appear to 
use the word "raw" to refer to the calibrated (filtered for 
3-year analyses, unfiltered for 5-year analyses) TOD and do 
not appear to have carried out the calibration step. 

Independently of these considerations, the correct way 
to compile the calibrated WMAP TOD into sky maps must 
clearly be the one that optimises the sharpness of the im- 
ages of bright sources. Whether or not the calibration step 
itself was carried out with a small timing error — leading 
to an artificial dipole-induced difference signal — remains to 
be investigated. Liu et al. (2010) 's calculations suggest that 
the latter is an interesting possibility. 
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Appendix A: Time offset between IVIeta Data 
Table and Science Data Table 

Since the time offset between the Meta and Science Data 
Tables in the WMAP TOD has escaped attention for many 
years, a method of verifying this is provided here. 

The TOD files, the GNU Data Language 
(GDL), and the IDL Astro library are available at 
http : / /lamib da . gsf c . nasa . gov/ product/map/ d.r2/, 

tod_f cal_ge t . cf m) 

http : / / idlastro . gsf c . nasa . gov/ ftp/ astron . dir . t ar . gzj 
http : //gnudatalcinguage . sourcef orge .net/, respec- 
tively. The file path for the IDL Astro library can be set 
up using the equivalent of 

cd pro && In -s */*.pro . && export GDL_PATH='pwd' 

in a bash shell. The line 

ndata = product ( axis, /integer ) 

in the IDL Astro routine fits/fits_READ.PRO should be 
replaced by 

ndata = floor ( product ( axis ) +0.5) 
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for use in GDL-0.9^rc3 or earlier. The following GDL script 
demonstrates the timing offset in the TOD, where "tod. fits" 
is the name of a TOD file. 

pro TOD_offset 

fits_open, 'tod. fits', handle 

; ; get table of times (metatime) of quaternions 
; ; of spacecraft attitude 
fits_read, handle, metatable, metaheader, $ 
extncune = 'Meta Data Table' 

tbinfo, metaheader, meta_h 

nietatinie= tbget(nieta_h, metatable, 'Time') 

; ; get table of times (scitime) stored in table 
; ; that includes observational data 
fits_read, handle , scitable , sciheader , $ 

extname = 'Science Data Table' 
tbinfo , sciheader , sci_h 
scitime= tbget(sci_h, scitable, 'Time') 

; ; re-interpret ID table as 2D table 

scitimeO=ref orm(scitime,30, 1875) 

; ; calculate the offsets in seconds 

offset = (metatime - scitimeO[0,*] ) *3600d *24d 

; ; print some examples 

print, 'quaternion time minus data time in s: '+ $ 

'first 10 frames of the day' 
print , offset [0:9], f ormat=" (dl5 . 7) " 
print, 'quaternion time minus data time in s: '+ $ 

'last 10 frames of the day' 
print , offset [1865 : 1874] , f ormat=" (dl5 . 7) " 
print, 'mean and stand, deviation:' 
print ,mean(of f set) ,f ormat=" (dl5 . 12) " 
print, sqrt(mean((off set-0.0256)~2)) , $ 

format="(dl5.12)" 

end 



